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Abstract. Motivated by Yang-Mills dark energy model, we propose a new model 
by introducing a logarithmic correction, we find that this model can avoid the 
coincidence problem naturally and gives an equation of state w smoothly crossing — 1 
if an interaction between dark energy and dark matter exists. It has a stable tracker 
solution as well. To confront with observations based on the combined data of SNIa, 
BAO, CMB and Hubble parameter, we obtain the best fit values of the parameters 
with lcr,2cr,3cr errors for the noncoupled model: = 0.276 ± 0.008tj]:J]}^t5^;5]^t 

h = 0.699 ± 0.003 ± 0.006 ± 0.008, and for the coupled model with a decaying rate 
7 = 0.2: rim = 0.291 ± 0.004tg °J]?t° gi2^ h = 0.701 ± 0.002 ± 0.005 ± 0.007. In 
particular, it is found that the non-coupled model has a dynamic evolution almost 
undistinguishable to ACDM at the late-time Universe. 
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1. Introduction 

The accelerating expansion has been supported by observations of supernova 
observations [HIS], cosmic microwave background radiation (CMB), and the large scale 
structure through the baryon acoustic oscillation (BAO). However, the physical origin 
of acceleration is still a challenging mystery. Within the framework of general relativity, 
this can be interpreted by a cosmic dark energy with negative pressure. The simplest 
dark energy model is the cosmological constant (ACDM), where the vacuum energy 
is responsible for the accelerating expansion. Whereas, it suffers from two problems 
degenerately. The first is the fine-tuning problem: The observed vacuum energy density 
of order ~ 10~^^ GeV^ is about 10^^^ orders of magnitude smaller than the value expected 
by quantum field theory for a cut-off scale being the Plank scale, and is still about 10^'' 
orders smaller even for a cut-off scale being the QCD scale [3]. The second is the 
coincidence problem: The conditions in the early Universe have to be set very carefully 
in order for the energy density of the vacuum and that of the matter to be comparable 
today. To solve these problems, abundances of dynamical dark energy models have 
been proposed, such as scalar field [H O El El [HI E], vector field [lOl [IH [121 HI], and 
decaying vacuum energy [HI |T5], Holographic dark energy model [El [TTl [13 [19] and 
so on (see [3] for details). In quintessence of scalar field, the Lagrangian density has 
a standard form £ = |(V0)^ — V{(f)), which contains a canonical kinetic term and a 
potential term. Some particular forms of the potentials of quintessence [H [5l [6] or the 
coupled quintessence with dark matter [201 [211 [22] , can not only lead to the late time 
acceleration of the Universe, but also avoid the coincidence problem. 

As well known, the quintessence always has an equation of state (EoS) w larger 
than —1. However, there have been some preliminary evidences that the current value 
of w may be less than —1, as indicated from observations [231 [211 [23 [211 [23 [23 [29] . 
On the other hand, the phantom field with a negative kinetic energy [30] predicts a 
EoS being always smaller than —1. Furthermore, another scalar field model with a 
noncanonical kinetic term called K-essence was proposed [HI [32l [33] . The Lagrangian 
is generally taken to be £ = /(0)(7(X), where X = |(V0)^. In [33], oc and 
g{X) is a polynomial of X, leading to w > —1 for a > as a quintessence and w < —1 
for a < as a phantom. For the considerations of the EoS of dark energy tracking that 
of the background, one expect that the EoS of dark energy may be lager than —1 in 
early times and smaller than —1 in late times, i.e., w crosses —1. If this is supported by 
further observations, many dark energy models would undergo a problem. Note that, 
even in K-essence model, w can not cross —1 for a fixed a. Is it possible to have a transit 
from quintessence to phantom? The particular interacting phantom dark energy could 
give smooth transit from w > —1 to w < —1 [M]. The non- minimal interaction between 
dark matter and dark energy with a single scalar field could also make possible to do 
the crossing of the phantom divide 05]. Besides, based on employment of two scalar 
fields, a class of models have been proposed, thereby using extra degrees of freedom 
[361 [371 [38]. 
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In this paper, we propose to study a new dark energy model described by a 
single scalar field, whose Lagrangian density contains a logarithmic factor. We call 
it effective scalar field (ESF) dark energy model. Since the kinetic energy term is 
noncanonical, it should belong to a subclass of K-essence models. This is inspired 
by our previous work on the quantum effective Yang-Mills condensate (YMC) dark 
energy model [IH |39l SQl SH ill |13] with C^ff oc F In F for 1-loop case [lOl [IH |39l HQ] , 
where F is the squared gauge field strength [3H HS], |T0] . The nonlinear kinetic terms 
appear generically in the effective action in string and supergravity theories [16]. The 
appearance of a logarithmic correction in the field is generic for effective quantum 
theories, e.g., the Coleman- Weinberg potential |17], the effective gravity ^Hl as 
well as the effective Yang-Mills field Since the nature of dark energy is still 

unknown, in the following, we will investigate the phenomenological properties of ESF 
model. Based on the observation of nearby galaxies [50], an interaction between dark 
energy and dark matter is favored since it would give a more rapid structure formation 
than predicted by the ACDM model [51]. So, we will also generally consider that an 
interaction between dark energy and matter exists. As will be seen, with one scalar field, 
the model provides a smooth dynamical transit from quintessence to phantom, with its 
w going from > — 1 at high reshifts to < —1 at low redshifts, if the field decays into 
matter. All the physical quantities involved in the model are smooth during the whole 
dynamical evolution. Moreover, the coincidence problem is also avoided in this model, 
since it has a stable attractor solution. But, unfortunately, the fine-tuning problem is 
still exist since a model parameter has to be tuned to accord with the low density of 
the dark energy component obtained through observations. We will demonstrate these 
afore-mentioned points, and also carry out a joint analysis for the model, confronting 
it with recent observations from SN la [53l HH |55] , BAO peak measurement of large 
scale structure from the Sloan Digital Sky Survey (SDSS) [56j and the Two Degree 
Field Galaxy Redshift Survey (2dFGRS) [57], the shift parameter of CMB [5H], and the 
history of the Hubble parameter [59l [60l [61] . Throughout this paper, we adopt a unit 
with c = 1. Greek indices /i, z/, ... range over to 3, and Latin indices i,j, ... range over 
1 to 3. 



2. The effective scalar field model 



We consider a spatially flat Universe described by Friedmann- Robert son- Walker metric 

ds"^ = dt^ - a'^{t)6ijdx'dx^ , (1) 

where the scalar factor a{t) is determined by the Friedmann equation: 
__n 8ttG , , 

= ^(P<^ + Pm + Pr), (2) 

where p^., pm and p^ represent energy density of dark energy, matter and radiation, 
respectively, and H = d/ai?, the Hubble parameter. The dark energy is described by a 
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scalar field with a Lagrangian density 

i(V0)2-V^(0) 



(3) 

ae 

where (V0)^ = g'^'^d^(j)d^(j), V{(j)) is a function of 0, a is a scale of energy density to be 
fixed by observations, and Ine = 1. in Eq.(^ is formally similar to the Lagrangian 
density of the 1-loop effective YMC dark energy model [11]. Assuming (p is homogeneous 
and isotropic, and only depends on time, i.e., = 0(t). The action is given by 

S = J d'xy/^C^, (4) 

where g = det{gf^^). The variation of the action (jl]) with respect to gives 

/a e\ ■ dV , , 

0+ 3- + - + — = 0, (5) 



a e 

where a dot denotes d/dt, and e = In (|0^ — V)/a. Note that, Eq.(|5]) differs from that 
of the quintessence model by the extra term {e/e)(f) . When a coupling exists between 
ESF and matter, their dynamical evolution equations are given by 

+ ?>H{p^ + p^) = -Tp^, (6) 

Pm + SHp^n = Tp^, (7) 

where F denotes the energy transformation rate from ESF to matter. For simplicity, 
we assume F is a constant. The radiation is an independent component and evolves as 
Pr(t) oc a~^. 

In the following we focus on the simple case of = 0. The energy density and 
pressure are easily gained by the variation of action (jlj) with respect to g'^'^: 

P0 = a{e + l)e^ p^ = a{e - l)e^ (8) 

with e = ln(i0^/a;) > —1 required by p<^ > 0. It is easy to prove that, for the non- 
coupled case (F = 0), Eq. reduces to Eq. with the help of Eq. ([H]). The EoS is 
given by 

w = — = -. (9) 

P^ e + 1 

In high energy limit with e ^ 1, w ^ 1, different from the high energy behavior w ^ 1/3 
of the YMC model [llj. At the critical point e = 0, one has w = —1. Furthermore, 
w < —1 will be arrived when — 1 < e < 0. Introducing dimensionless x = Pm/d and 
r = Pr/ct, Eqs. (E]) and (I7j) read as 

.' + -^ + 4^ = 0, (10) 
e + 2 g{e + 2) ' ^ ^ 

x' + 3x--{l + e)e' = 0, (11) 
Q 

where ' = d/dN with = \na{t), 7 = F/ (^ ^^^" ^ ^ is the dimensionless decaying rate, 
and g = [x + r + e^(e + 1)]^''^. Given initial values (£j, Xj), Eqs. ( TTUj) and (ITT]) can be 
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solved for each 7. 7 > means that dark energy decays into matter, and vice versa. 
In this paper, we assume that 7 is positive, since 7 < will lead to a negative matter 
density in the future, which is unacceptable from the point of view of physics. 

Firstly, we would like to discuss the simple case of the non-coupled (7 = 0) ESF 
model. To ensure the standard cosmology not being spoiled by the presence of dark 
energy, one may take p^ij Pri < 10~^ at Zi ~ 10^. The outcome is that, for any Si in 
a wide range — 1 < < 55 corresponding to p^i ranging over almost infinity orders 
of magnitude, the current status (fi<^, — (0.73, 0.27) is always attained. So the 
coincidence problem is solved at the price of choosing a fixed a. As can be seen in FigJT] 
(a), different initial values of Si = —0.9,0, 10,50, lead to the same density of the ESF 
at present time. During earlier stages the decreasing P(j,{z) is subdominant to Pr{z) and 
Pm{z)- Note that P(i,{z) always levels off at a certain time, earlier for a smaller Si. Then, 
it surpasses pr aX. z 10 and surpasses p™. at 2; ~ 1, respectively. Whereas, the matter 
density Pmif) evolves independently as Pmii) oc since it does not couple with the 

ESF. FigH] (b) shows that w will decrease and increase with time and approach —1 at 
the late-time for £j > and Si < 0, respectively. Moreover, for Ei = 0, the ESF acts as 
the ACDM model exactly. Note that, in the above three cases, all the corresponding w 
will stay at —1 in the future and never cross —1. It can be understood as follows. For 
7 = 0, Eq. ( !T0|) has a solution: 

= Aa-\ (12) 

where A is an integration constant determined by the initial conditions. Except the 
particular case that e = 0, A is always positive, i.e., the RHS of Eq. ( fT2|) will always 
be larger than 0. This indicates that e will never cross no matter the initial value 
of e is positive or negative. According to Eq. (Q, we know that w can not cross —1 
resulting from e failing to cross 0. In the special case of Ei = 0, A is fixed to be zero. 
Thus, e = will be kept all the time, i.e., w = —1 is constant. Since e'^e^ decays with 
the expansion Universe as oc a~^, e — )• at ^ = 0. Therefore, p^{z = 0) ~ a ~ 0.73pc, 
where pc is the critical density of the Universe. We find that, for the whole range of e^, 
the resulting dynamical evolution in the recent past {z < 30) is almost identical to that 
in ACDM with deviations (p^ — pa)/pa < 10^^. The total EoS is Wtot = J2 ^j'^j^ where 
j stands for the ESF, matter and radiation, respectively. In the future (A^ — ?■ 00), one 
has Qr ^ 0, Qm — ^ and ^ 1. Hence, Wtot = w = —1, as shown in FigH] (b). That 
is, the Universe will do an exact de Sitter expansion, and there is no big rip event which 
some dark energy models would encounter. Note that, the interacting phantom dark 
energy could also avoid the big rip event [M] . 

Secondly, as an exphcit example of ESF, we discuss the coupled case with 7 = 0.2. 
The initial conditions for the ESF and radiation are chosen the same as the case of 
non-coupled ESF. The initial condition for the matter is chosen a little differently from 
the non-coupled case in order to ensure the current status {Q(f,,Qrn) — (0.73,0.27). 
As illustrated in Figj2] (a), the coincidence problem is also avoided in this case. The 
dynamic evolutions of the density of the coupled ESF are quite similar to those in the 
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Figure 1. (a) p^{z), pm{z), and Pr{z) in tlie non-coupled ESF model. For various 
initial values Si the current status VL^ = 0.73 and VLm = 0.27 is always achieved, (b) 
w{z) does not cross —1 without the coupling. 



non-coupled ESF. However, levels off around 2; ~ and will approach a constant 
instead of decaying as oc a{t)~^ . This is caused by the coupling 7 7^ 0. Fig. Et^b) plots 
evolutions of the corresponding w{z). Due to the coupling, w crosses —1, arrives at 
wq ~ —1.05 at present, and settles down to a constant value w ~ —1.07 in future. The 
influence of the coupling 7 has been investigated, and computations show that a greater 
7 yields a larger matter fraction Vt^ and a smaller EoS Wq at present. We have also 
found an interesting relation: VL^ = —1/w a,s N 00, which is similar to YMC model 
[62] . This implies that the total EoS satisfies Wtot = ^,f,w = — 1 as — )■ 00. Thus, the 
coupled ESF model also predicts an exact de Sitter expansion in future, and the big 
rip event is avoided naturally. The parameter a in this case are determined by setting 
P0 ^ 0.73pc at z = 0, leading to a ~ 0.76pc. Unfortunately, the particular choices 
of a being the same order of magnitude as pc let the ESF model still suffer from the 
fine-tuning problem. 

We have carried out an analysis of dynamic stability of the set of Eqs. ( fTOll 
and (fTTj) . and found that it has the fixed point {Sc^Xc) = (0,0) for 7 = and 
{EcXc) = (—0.0323,0.0625) for 7 = 0.2, respectively, as — )■ cxd. Moreover, any 
perturbations 6e and 6x are both decay as e~^^ for 7 = 0, and decay as linear 
combinations of e~^-'^^^^^ and e"^-^^''^^ for 7 = 0.2, respectively. Thus, the fixed points 
are stable, and then the attractor solutions of p</,(t) are obtained for the above two cases. 
Aside the trajectory smoothness in phase space, the stability is problematic in single 
scalar field models [63] • Moreover, all the physical quantities in our models, such as 
a(t), P(f,{t), and w{t), are smooth from the initial moment up to the future. 
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Figure 2. (a) p^{z), Pm{z), and Pr{z) in the coupled ESF model with 7 = 0.2. 
For various initial values the current status = 0.73 and fi™ = 0.27 is always 
achieved, (b) w{z) in the coupled model. Due to coupling, w{z) is able to cross —1 at 
low redshifts. 



3. Constraints from SN la, BAO, CMB and Hubble data 

For a model of dark energy to be viable, it needs to confront or be constrained with 
observational data. Here we constrain the V = model with the latest observational 
data of the 557 SN la assembled in the Union2 compilation [55] , the BAO measurement 
from SDSS [56] as well as 2dFGRS [57], the shift parameter of CMB from WMAP7 [58], 
and the the history of the Hubble parameter [591 EOl l6T| . 

First, we compare the theoretical distance modulus to the observed ones compiled 
in [55]. The theoretical distance modulus is defined as 

fithiz) = 5 logio + /io, (13) 

where Dl{z) = Hq{1 + z) Jq is the Hubble-free luminosity distance in a spatially 
fiat Universe, and /io = 42.38 — 51og]^o^ with h the Hubble constant in the unit of 
100km/ sec/ Mpc. The late-time Hubble rate of the effective scalar model is given by 

H{z)=Ho[Q^{z) + Qm{^)f', (14) 

where Qi{z) = pi{z)/pc for i = (j),m. Since the evolution for z < 10^ in this model is 
insensitive to the initial conditions, we choose Si = 10 in the following calculations for 
concreteness. For the SN la data, the function is 

x|.(P.;/^o) = f: ^^"^"'^~^"^"-^^' , (15) 

1=1 
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where Ps stands for a set of parameters, such as Vtm- The nuisance parameter /iq can 
be analytically marginalized over [61], so that one actually minimizes xIatIPs) instead 

of x|jv(Ps;/io)- 

Next, the BAO is revealed by a distinct peak in the large scale correlation function 
measured from the luminous red galaxies sample of the SDSS a.t z = 0.35 |56], as well 
as in the 2dFGRS at ;z = 0.2 [57]. The peaks can be associated to expanding spherical 
waves of baryonic perturbations. Each peak introduces a characteristic distance scale 



Dv{zbao) 



zbao 
H{zbao 




1/3 



(16) 



The observational date from SDSS and 2dFGRS measurements yield D^(0.35) / D^{Q.2) = 
1.736 ± 0.065 |[57j. The best fit values for the model are given by minimizing [SH E2] 

2 . , _ {[D,i0.35)/D,{0.2)]th - [I^.(0.35)/D,(0.2)]ofe,)^ 

XbAo[Ps) — —2 ; U'J 

'-' Dy (0.35)/ Dy (0.2) 

where a|5^(o.35)/D4o.2) = 0-065. 

As discussed in [661 EI], the first peak of the CMB spectrum of anisotropies, h, is 
more suitable to be used to test the interacting dark energy model than the CMB shift 
parameter, R = ^/^l^D L^Zrec) / + Zrec) [HHl [69], where Zred = 1091 [58] is redshift of 
recombination. Then, we use /i, which is related to the angular scale, I a, by [70] 



^1 = - Si] 



where 



p ^0.1 



0.267 ^^j (19) 

with p = Pr{zrec) / Pm{zrec) the density ratio of radiation and matter at the time of 
recombination. The acoustic scale is defined as 

= ^ / TTTt/ / ITTT' (20) 

Jo n [Z) J Zrec rl [Z) 

where the sound velocity is = (3 + ^^), with VL^ and VL^ the present density 
parameters of baryons and photons, respectively. With the observed position of the 
first peak hobs = 220.8 ± 0.7 [H], the for CMB is 

x'cMBiPs) = i^H^^i^^, (21) 
where o"; = 0.7. 

Finally, the Hubble parameter as a function of redshift z can be written as 

Hiz) = (22) 

^ ' l + zdt ^ ^ 

Then, once dz/dt is known, H{z) is obtained directly. Simon et al. [72] and Stern et al. 
[59] obtained H{z) in the range of < 2; < 1.8, using the differential ages of passively- 
evolving galaxies and archival data. Recently, some high precision measurements 
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Figure 3. Left: The confidence contours for the pair of free parameters (r2m,/i) 
obtained by constraining the non-coupled ESF model with the joint observational data 
from SN la, BAO, CMB and H(z). Right: The confidence contours of the same pair of 
parameters of the coupled ESF model with 7 = 0.2 obtained by the same observational 
data. 



constrained H{z) aX z = from the observation of 240 Cepheid variables of rather similar 
periods and metallicities ^U\. Besides, H{z) at z = 0.24, 0.34 and 0.43 is obtained [HI] 
by using the BAO peak position as a standard ruler in the radial direction. We employ 
the twelve data in [60l [59] and the three data in [61]. The best fit values of the model 
parameters from observational Hubble data are determined by minimizing 

xL>(P.)-i: "'"-'^-' (23) 

i=l ^ K^i) 

Thus, the total is combined as 

Xtotal — XsN + XbAO + XcMB + Xllub- (24) 

As the likelihood function is determined as £ oc exp{—Xtotai/'^)y the best fit values of 
flm and h follow from minimizing Eq. fl24|) . Fig|3] shows the 68.3% (Icr), 95.4% {2a) and 
99.7% (3cr) confidence contours in the Qm — h plane for both the ESF model with 7 = 
and 7 = 0.2. For 7 = 0, the best fit values of 1-dimension up to 3a confidence level 
are: fi^ = 0.276 ± 0.008l|]:g}^;[|;[j^^, h = 0.699 ± 0.003 ± 0.006 ± 0.008, with a minimal 
Xl^^^i = 554.713; while for 7 = 0.2, the results are: = 0.291 ± 0.004t[!:[]J]^l[]:{]}?, 
h = 0.701 ± 0.002 ± 0.005 ± 0.007 with Xtotoai = 556.033. For comparison, we also 
calculate the case of ACDM and find it gives almost the same results as the non-coupled 
ESF model at a very high precision. 

Now, we would like to compare these three dark energy models, i.e., ACDM, ESF 
models with 7 = and 7 = 0.2. A conventional criterion for comparison is Xmin/dof, in 
which the degree of freedom dof = N—k, whereas N and k are the number of data points 
and the number of free model parameters, respectively. We calculated the xlninl^^f for 
the three models, which can be seen in Table. I. Besides, there are other criterions 
for model comparison such as the Bayesian evidence [731 HI]- However, the Bayesian 
evidence is usually sophisticated. As an alternative, we can use some approximations 
of Bayesian evidence such as the so-called Bayesian Information Criterion (BIC) and 
Akaike Information Criterion (AIC), instead [75]. The BIC is defined as [76] 

BIC = -21n£^a^ + fclnA^, (25) 
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Table 1. Comparison of the three models considered in this work. 



Model 




^total 


xVdof 


ABIC 


AAIC 


ACDM 




554.713 


0.970 








ESF (7 


= 0) 


554.713 


0.970 








ESF (7 


= 0.2) 


556.033 


0.972 


1.32 


1.32 



and AIC is defined as [77] 

AIC = -21n£^,, + 2fc, (26) 

where C^ax is the maximum hkehhood. In the Gaussian cases, xLin = ~21n£max- So, 
the differences of BIG and AIG between two models are ABIG = Axmin + ^ 
and AAIG = Ax^m + 2AA;, respectively. In Table. I, we also present the ABIG and 
A AIG. One can find easily from Table. I that, the non-coupled ESF model and AGDM 
model not only have an almost identical evolution in the recent past {z < 30), but 
also are undistinguishable in confronting with the combining observations from SN la, 
BAO, GMB and Hubble parameter. Moreover, for the coupled ESF model, it only gives 
a little larger Xmin little larger values of all the criterions for model comparison. 
Thus, 7 > would be favored if further observations support w < — 1 as indicated 
in Refs. [25], since the coupled and non-coupled ESF models perform similarly in 
analysis. 

4. Conclusions 

Inspired by a generic feature of effective quantum fields, we have proposed a scalar 
field dark energy model, whose Lagrangian contains a logarithmic correction. It can 
be regarded as a special case of the generic K-essence models. For an initial value p^j 
ranging over almost infinite orders of magnitude, p^{t) tracks radiation and matter, 
and the current status ^1^) — (0.73, 0.27) is always attained. So the coincidence 
problem is solved if the parameter a is chosen in advance, but the fine-tuning problem 
remains. Moreover, w smoothly crosses —1 if the ESF decays into matter. For a decay 
rate 7 = 0.2, the EoS arrives at wq ^ —1.05 at present. A greater 7 yields a larger Qm 
and a smaller wq at present. As t — ?■ 00, the expanding spacetime approaches the de 
Sitter as an asymptote, which is also a stable attractor, and there is no cosmic big rip. 
For the non-coupled model, w approaches —1 but does not cross —1, and the dynamic 
behavior is almost the same as AGDM for low redshifts. In particular, for an initial 
Ei = 0, the model reduces to AGDM. Since the meaning of a non-zero V{(j)) is unknown, 
we did not discuss the properties of the ESF model with V^cj)) 7^ 0. Some particular 
forms of V{(j)) would be investigated in the future study. 

In confronting with observations of SN la, BAO, GMB and Hubble parameter, 
we plotted the confidence contours in the — h plane for the ESF model with 
7 = and 7 = 0.2. The best fits of the parameters are: flm = 0.276 ± 0.008 
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and h = 0.699 ± 0.003 with xLoai = 554.713 for 7 = 0; fi^ = 0.291 ± 0.004 and 
h = 0.701 ± 0.002 with Xtotoai ~ 556.033 for 7 = 0.2. Furthermore, the non-coupled 
ESF model is distinguishable from ACDM model under present observations. Besides, 
we compared the three dark energy models studied in this work using x^/dof, BIG and 
AIC . It is found that a non-coupled ESF model is a little more favored, however, the 
coupled model will survive if further observations support w < —1 strongly. 
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